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ABSTRACT 

We examine the proposal that the subset of neutron-star and black-hole X-ray binaries 
that form with Ap or Bp star companions will experience systemic angular-momentum 
losses due to magnetic braking, not otherwise operative with intermediate-mass com- 
panion stars. We suggest that for donor stars possessing the anomalously high mag- 
netic fields associated with Ap and Bp stars, a magnetically coupled, irradiation-driven 
stellar wind can lead to substantial systemic loss of angular-momentum. Hence these 
systems, which would otherwise not be expected to experience 'magnetic braking', 
evolve to shorter orbital periods during mass transfer. In this paper we detail how 
such a magnetic braking scenario operates. We apply it to a specific astrophysics 
problem involving the formation of compact black-hole binaries with low-mass donor 
stars. At present, it is not understood how these systems form, given that low-mass 
companion stars are not likely to provide sufficient gravitational potential to unbind 
the envelope of the massive progenitor of the black hole during a prior "common- 
envelope" phase. On the other hand, intermediate-mass companions, such as Ap and 
Bp stars, could more readily eject the common envelope. However, in the absence of 
magnetic braking, such systems tend to evolve to long orbital periods. We show that, 
with the proposed magnetic braking properties afforded by Ap and Bp companions, 
such a scenario can lead to the formation of compact black-hole binaries with orbital 
periods, donor masses, lifetimes, and production rates that are in accord with the 
observations. In spite of these successes, our models reveal a significant discrepancy 
between the calculated effective temperatures and the observed spectral types of the 
donor stars. Finally, we show that this temperature discrepancy would still exist for 
other scenarios invoking initially intermediate-mass donor stars, and this presents a 
substantial unresolved mystery. 
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1 INTRODUCTION 

Magnetic braking is a relatively poorly understood mech- 
anism for angular momentum loss in single stars as 
well as compact binaries, but is often invoked as the 
driver of mass transfer in compact binaries such as cat- 
aclysmic varia bles and l ow-mass X-ray binaries (see, e.g., 
Verbunt. fc Zwaant.l98ll : iRappaport. Verbunt fc JossHlQSa 
Spruit fc Rittejll983^ . Models of magnetic braking usually 
invoke two key components: a large-scale magnetic field and 
a stellar wind. When coupled together, the field provides a 
long lever- arm which makes the wind an efficient angular- 
momentum sink. 

According to conventional thinking, braking via a mag- 
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netically constrained stellar wind is inoperativ e in stars 
with radiative envelopes, (see, e.g. lKawaleJll988h and thus 
is confined to stars with M < 1.5 Mq. However, a sub- 
class of intermediate-mass stars (the Ap and Bp stars) pos- 
sess substantial magnetic fields; herein we suggest that an 
irradiation-driven stellar wind combined with their magnetic 
fields would lead to significant magnetic braking in compact 
binaries containing such stars. 

We apply this heretofore unutilized magnetic braking to 
the formation of compact black-hole X-ray binaries. These 
systems are not easy to account for within previous forma- 
tion scenarios - yet they comprise about half of the well- 
studied black- hole X-ray binaries. Nine of th e seventeen sys- 
tems listed in iLee. Brown fc WiiersI ^|2002^. iRitter fc Kolbl 
(:2003) and iPodsiadlowski. Rappaport fc Hani tOOA) (here- 
inafter PRH, and references therein) have orbital peri- 
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ods < Id and inferred dono r masses < 1 (see also 
McClintock fc Re^nillar'dll2005^■ The empirical estimates of 
WiiersI l)l99fih and lR,omanil tmff) suggest that the Galactic 
population of short-period black-hole binaries may number 
over 1000 systems. 

The standard formation scenario for such compact 
black-hole binaries typically requires a low-mass donor to 
have ejected the envelope of its primordial binary com- 
panion - the massive black-hole progenitor - as it spi- 
rals inwards during a common-envelop e phase (see, e.g. 
Ide Kool. van den Heuvel fc Pvlvseill98'il also PRH and ref- 
erences therein). The orbital period is reduced from years to 
days, and the core of the massive star continues to evolve 
toward collapse and the formation of the black hole in 
the system. However, it is not clear how the energy re- 
quired to unbind the envelope of the black-hole progeni- 
tor can be provided by a secondary star with a mass of 
only < 1.5 Mp, (e.g. [Podsiadlowsk i. Cannon fc Reejll995l: 
iPorteeies Zwart. Verbunt fc Ergmal Il997l: iKaloeeral ll99S . 
PRH; see the Appendix for a brief summary of this argu- 
ment). In the conventional scenario, this problem cannot 
simply be resolved by assuming a more massive secondary. 
Donor stars with M > 1.5 Mq are not expected to evolve 
into the short-period population, as they are thought not 
to be subject to canonical magnetic braking; mass transfer 
from the less massive to the more massive component of the 
system tends to widen the orbit (a consequence of conserva- 
tion of angular- momentum) . 

In addition to this energetics problem, standard as- 
sumptions about star formation also suggest a demographic 
difficulty - even if the envelope of the massive primary could 
be ejected by the spiral-in of a low-mass star, there are ex- 
pected to be relatively few primordial binaries with such an 
extreme mass ratio as required in that standard pict ure (e.g. 
[Garmanv, Conti & Masscv 1980; Fisher, Schroder fc SmithI 

Hence, attempts to produce these low-mass black- 
hole binary systems have sometimes been exotic, pos- 
tulating _^^oi_exam2lB-_descent from a triple sys- 
tem (Eggleton fc Verbunt|[l98§) . that the companion and 
black hole both form from a Th orne-Zytkow object 
jPodsiadlowski. Cannon fc Reesll995ll , or as a result of run- 
way accretion in the merger of a neutron-star with a massive 
star (PRH) . Such ideas have been considered since the 'sim- 
ple' alternatives often have far-reaching consequences, e.g., 
modifying the energetics of common-envelope evolution by 
a large factor to allow a low-mass star to eject the envelope 
of the primary is unpalatable, whilst magnetic braking from 
intermediate-mass radiative stars is not otherwise expected. 

For standard input physics, PRH concluded that 
"black-hole binaries with low-mass secondaries can only 
form with apparently unrealistic assumptions", which is 
an unsatisfactory position for binary evolutionary theory. 
PRH demonstrated that for the expected range of common- 
envelope parameters, black-hole binaries are much more 
likely to form with intermediate-mass companions than low- 
mass companions. They also pointed out that to form the 
observed short-period population from these intermediate- 
mass systems, an additional angular-momentum loss mech- 
anism, such as anomalous magnetic braking, would be re- 
quired. 

An elegant way to solve this problem of forming short- 



period black-hole binaries would be to find such an angular- 
momentum loss mechanism for intermediate-mass stars. Ide- 
ally it would operate only on a subset of the intermediate- 
mass donor population, such that the long-period and short- 
period black hole binary populations can both be formed. 
This paper suggests a self-consistent and plausible way for 
initially intermediate-mass companions to lead to the low- 
mass companions currently observed, at both long and short 
orbital periods. 

Though the bulk of intermediate-mass stars are not ex- 
pected to undergo magnetic braking, a fraction of such stars 
are known to have anoma lously stronp magnetic fields, the 
so-called Ap and Bp stars IMoss 1989: ' Braithwaite fc Spruil 
12004") with surface field strengths ranging from 100 G to over 
10,000 G. When such a magnetic field is combined with an 
irradiation-driven wind from the surface of the donor star, 
a plausible angular-momentum loss mechanism is obtained. 
This paper is motivated by the realisation that such a mech- 
anism can act on a portion of the intermediate-mass pro- 
genitor population to produce the short-period low-mass bi- 
nary companions to black holes. Only a small portion of the 
intermediate-mass donors are affected - those with anoma- 
lously high field strengths. As we will show in Section 3, 
the lifetimes of systems with high magnetic fields are sig- 
nificantly longer, since the initial high mass-transfer rate 
quickly reduces the mass of the donor stars, making them 
evolve as lower-mass stars with much longer evolutionary 
timescales. This increases the relative fraction of the short- 
period black-hole binary population. 

We find that the proposed magnetic braking mechanism 
should function as anticipated, driving systems to short or- 
bital periods and producing low-mass donors. However, we 
also find that low-mass donor stars which descended from 
intermediate-mass progenitors have rather higher temper- 
atures than can be consistent with the spectral types ob- 
served for the donor stars in short-period black-hole X-ray 
binaries. If no reconciliation of this discrepancy is found, 
this result is itself potentially very illuminating, as it would 
apply to any formation channel which invokes primordially 
intermediate-mass donor stars. This would imply that ei- 
ther our understanding of common-envelope ejection needs 
severe revision, or an "exotic" channel is responsible for the 
formation of short-period black-hole binaries. 

In section 2 we assemble the relevant physics to quan- 
tify the proposed anomalous magnetic braking mechanism 
and investigate its properties. In Section 3 we demonstrate 
how this magnetic braking naturally produces the observed 
short-period black-hole binaries. Finally, we consider obser- 
vational tests of the model and illustrate its successes and 
major shortcoming. 



2 ASSUMPTIONS AND DERIVATIONS 

Magnetic braking removes rotational angular-momentum 
from a star when its stellar wind remains co upled to the stel- 
lar magnetic field as it leaves the star (e.g. I Weber fc Davis! 
Il967t iMestel fc Spniilll987l: iKawaleJ Il98a) . We make the 
standard assumption that the wind corotates out to the 
magnetospheric radius, defined as the radius where the mag- 
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(1) 



netic pressure and ram pressure are balanced, i.e., 

1 ^2^B^^ B^B^ 

2^" ~ Stt Stt r6 ' 
where r is the radial distance from the centre of the star, 
p and V are the density and velocity of the stellar wind, 
respectively, Bs is the magnetic field strength at the surface 
of the star, and Rd is the radius of the donor star. Here 
we have assumed a dipolar magnetic field structure. We can 
then combine eq. (1) with the continuity of mass applied to 
the stellar wind- loss rate, Mwind: 

2 



(2) 



and the assumption that the wind velocity at the surface of 
the star is of order the escape speed: 



2GMd 



Rd 
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(where Md is the mass of the donor star under consideration) 
to yield the magnetospheric radius: 
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We note that assigning a velocity to the wind by assuming 
it is equal to the escape speed at the magnetosphere, rather 
th an the surface of the star, re covers the expression due 
to iLamb. Pethick fc Pine3 il973l) . We have performed the 
derivations below using both versions of the magnetospheric 
radius, and the differences in outcome are minimal. 

The rate of change of angular-momentum due to mag- 
netic braking, Jmb, is then: 



Jmb 



-1/4 



(5) 

where Qd is the angular rotation frequency of the donor star, 
and we have substituted in the expression for from eq. 
(4). We expect that Qd ~ ^oih for close binaries as they are 
likely to be tidally locked. 



2.1 Estimate of the Required Wind-Loss Rate 

Before proceeding, we first estimate what angular- 
momentum loss rate, Jmb, will prevent the mass transfer in 
the binary (from the lower-mass donor to the higher-mass 
black hole) from widening the orbit and, in fact, will allow 
for orbital shrinkage. To make a rough estimate of the effect 
of Jmb on the orbital separation, we assume conservative 
mass transfer and, from angular-momentum conservation, 
obtain: 



Jmb 



a Mrlof I ^ 
a Md 



Md \ 
Mbh / 



(6) 



where Mrlof is the mass-loss rate from the donor star that 
is transferred to the black hole via Roche-lobe overfiow, and 
a is the orbital separation. If we require that a < 0, and solve 
for Jmb , neglecting the mass of the donor star compared to 
that of the black hole, we find: 

Jmb > Mrlof (GMtotfl)''" , (7) 

where Mtot is the total mass of the binary. When we combine 
this expression with the middle term in eq. (5) and assume 
that the donor star is corotating with the orbit, we find: 
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If we now parameterise the orbital separation in terms 
of the donor st ar radius and t he dimensionless Roche-lobe 
radius, (e.g. lEggletonlirQSSTl as a = Rs/tl, and use the 
magnetospheric radius from eq. ^ we find the requirement 
for Mwind: 
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For typical mass ratios of interest, rz, — 1/3, and a donor 
star mass-radius relation of approximately Md oc R^^^^ (for 
CNO cycle homology), we have: 



/ Mwind 



V M© yr- 



> 5 X 10' 



_}_( Md\ Mrlof 



V Mq yr- 



(10) 

where Bs is expressed in Gauss. 

This yields an estimate of the required mass-loss rate 
in the magnetically coupled wind of Mwind ~ 4 X 10-'° Mq 
yr~^ for the illustrative parameters: Bs = 1000 G, Md = 
5 Mq, and Mrlof = 10"* Mq yr"'. Though our compu- 
tational work does not explicitly employ eq. (10), it helps 
illustrate the magnitude of an irradiation-driven wind that 
is required to significantly affect the evolution of binary sys- 
tems. 



2.2 Irradiation-Driven Winds 

It has been suggested that a substantial stellar wind may be 
driven from the donor star in a compact binary by the flux 
of X-radiati on that is produced by accretion onto the co l- 
lapsed star iRuderman et al ]|l98gl : lTavani fc LondoJl993ll . 
Using an even smaller wind-driving efficiency than suggested 
by Tavani & London, we will show that anomalously mag- 
netic intermediate-mass donor stars in close binaries can be 
braked sufficiently to form short-period systems. We will 
also show that for the majority of donor stars, i.e., ones 
with relatively weak magnetic fields, magnetic braking via 
an irradiation-driven stellar wind is not strong enough to 
drive the systems to shorter periods. In this same con- 
text, we note that llben. Tutukoy fc YungelsonI il995ll and 
llben. Tutukov fc Fedoroval Jl997l) have also suggested that 
an induced stellar wind could constitute an important dif- 
ference between the evolution of cataclysmic variables and 
low-mass X-ray binaries. Iben, Tutukov fc Fedorova further 
pointed out that a radiation-enhanced stellar wind should 
have an effect on the magnetic braking, but their calcula- 
tions did not examine the consequence of this, and they con- 
fined themselves to the orbital angular-momentum carried 
away directly by the wind - i.e. the specific orbital angular- 
momentum of the donor star. 

The stellar wind-loss rate is obtained by assuming that 
a fraction of the accretion luminosity is converted into the 
kinetic energy of a wind, such that the matter in the wind 
becomes marginally unbound (consistent with our assump- 
tions about the wind velocity). This energy-balance argu- 
ment gives: 



Mw 



Lxfnft 



Rd 
GMd 



(11) 



where the total accretion luminosity Lx is multiplied by a 
geometric factor fa to find the flux that intercepts the donor, 
and a wind-driving energy efficiency factor, f^. The X-ray 
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luminosity is taken as some fraction, /a; , of power produced 
by an accretion rate of Mrlof before the last stable orbit 
around a non-rotating black hole: 



V8 
3 



Mrlof c 



(12) 



where c is the speed of light, fx is a rest-mass to energy 
conversion factor of order unity, and the factor in brackets is 
the dimensionless energy lost by the innermost stable orbit, 
i.e., at three Schwarzschild radii. Finally, we find: 



RdMuhOF 



GMa 

where, for convenience, we have defined: 



(13) 



(14) 



Initially we used a combined value for tp/c^ of ~ 10"'', 
based on the assumptions that the wind-driving energy con- 
version efficiency /e ~ 10""^, the solid angle /n — 10~^ and 
/x ~ 1. Clearly all these individual values are subject to 
some uncertainty, but we believe t hat the assumed values o f 
tp are on the conservative side; e.g. lTavani &: LondonI lll993ll 
calculated 'wind efficiencies' - ft - between 10~^ and 10~^. 
We began by using their lowest calculated values for , and 
then took tp to be one and two orders of magnitude lower 
(i.e., ip/c^ = 10"^ and 10"^). 

In the computational models we took into account the 
fact that the irradiation-induced stellar wind carries away 
the specific orbital angular-momentum of the donor star in 
addition to the braking effect it produces on the rotation of 
the donor. We find that the direct loss of orbital angular- 
momentum via the stellar wind is quite small, as might be 
expected from the fact that the Alfven radius is much greater 
than the dimensions of the binary system. 
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Figure 1. The equilibrium mass-transfer rate solution predicted 
by eq. 1181 for a constant accretor mass of 7 Mq. We plot three 
different stellar magnetic field strengths - 300 G (dashed), 1000 G 
(unbroken line), and 3000 G (dotted). Here i/i/c^ = 10"*^. 
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2.3 Analytic Results for Mrlof 

A more sophisticated estimate than eq. llOt for the 
irradiation-induced wind loss can now be assembled. In par- 
ticular, we can form a closed set of equations, eliminate 
Mwind, and compute the mass-transfer rate when the mag- 
netic braking torque just prevents the orbital separation 
from increasing due to the effects of mass transfer. Begin- 
ning with eq. I^l, we substitute for the wind- loss rate from 
eq. (13), and use Kepler's 3rd law to obtain: 



Jmb ~ —Be 



V'MRLOFMtot / Rf \ 

a3 GM| ) 



1/4 



(15) 



The time derivative of the expression for the total 
angular-momentum of the binary. 



Jsys = MdMe 



Ga 



Mu 



(16) 



where Mbh is the black hole mass, with the assumption of 
mass conservation (i.e., Mtot = 0) and the constraint that 
d = (i.e., we are seeking the condition where the binary 
becomes neither wider nor tighter) yields: 



J sys — 



Ga 



Mu 



(Md - Mbh) Mrlof. 



(17) 



Figure 2. We show the evolution of two LMXB systems with 
initial donor masses of both 1 and 1.5 Mq, where each system 
contains a 1.4 Mq NS accretor. We plot the donor mass against 
time; the initial orbital periods were 8 h and 12 h, respectively. 
The solid lines use a standard magnetic braking parameterisation, 
whilst the broken lines use our new formalism with a constant field 
strength of 300 G (dot-dashed) and 500 G (dashed). 



Thence we can equate eqs. I15II and 11711 and solve for 
Mrlof. We again use the parameterisation that a = Rs/rL 
and define q — Md/MsH to find: 



Mrlof = i'S'^rl 



GHIi 



q + 1 

Q-1 



(18) 



FigureQillustrates the behaviour of this equation, again 
assuming Rd oc . Note that the solutions scale linearly 
with the overall wind-driving efficiency tp, and go as the 
square of Bs- 



2.4 Effects on the Canonical LMXB Population 

Before applying this model to Ap-star donors - adding 
an irradiation-driven wind to their intrinsically high field 
strengths - we did a calculation to check that such 
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physics would not unreasonably distort the standard pic- 
ture of low-mass X-ray binary (LMXB) evolution dur- 
ing the X-ray phase. These systems contain neutron-star 
accretors and low-mass donor stars, and they are typi- 
cally assumed to evolve via ' c onventional' magnetic brak- 
ing JVerbunt. fc ZwaanI Il98ll iRappaport. Verbunt fc Josj 
119831 ) where the magnetic field of the donor is presum- 
ably generated by dynamo action. Stellar dynamos are 
by no means precisely-so l ved p henomena, but following 
ICollier Cameron fc Jiankg lll994) we assume that the dy- 
namo of the donor star is saturated at the orbital periods 
we consider 1 d). Hence we use fixed stellar dipole field 
strengths of 300 and 500 G. 

For this simple assumption we found a range of LMXB 
evolutionary sequences that qualitatively reproduced test 
calculations performed with the magnetic braking parame- 
terisation of[RaDDaDort, Verbunt & ,]oss (1983). For clarity, 
the systems subject to our new mechanism do not experi- 
ence conventional magnetic braking, though we do include 
gravitational wave radiation. Figure |5| contrasts two of these 
binary sequences: the irradiation-driven braking is, for these 
assumptions, less strong than conventional magnetic brak- 
ing, and - when combined with conventional magnetic brak- 
ing - does not greatly affect the standard picture of LMXB 
evolution. Note that in this case we assumed tp/c? = 10"'' - 
the largest value of our wind-driving parameter. 



3 BLACK-HOLE BINARY POPULATIONS 

For populations of both short-period and long-period black- 
hole binaries to be produced, their progenitors must some- 
how differ. Something must distinguish between those sys- 
tems which are to evolve to long and to short periods. Many 
of the donor stars in black-hole binary systems are initially 
expected to be A and B stars (see, e.g., PRE). We there- 
fore suggest that a natural distinction, within a population 
of intermediate-mass stars, is between those which possess 
anomalously high magnetic field strengths (i.e., the Ap and 
Bp stars) and those which do not (ordinary A/B stars). The 
magnetic braking we postulate, that takes place in the high- 
field systems via irradiation induced stellar winds, would 
cause a subset of the black-hole binaries to evolve to short 
orbital periods. 

The question of what happens to Ap star field strengths 
when Roche-lobe overflow begins is not clear, and in this 
work we assume that the initially strong field persists. If 
the flux is a fossil field fro zen into the stellar matter (e.g. 
iBraithwaite fc Spruitl2004l) . the anomalous field may be lost 
with the mass in the envelope as it is transferred to the 
black hole. Note, however, that the field has to persist only 
long enough for a dynamo to begin operating as the stel- 
lar mass decreases. There has al s o bee n a suggestion by 
iTout. Wickramasinghe fc Ferrariol i2004ft that the magnetic 
fields seen in white dwarfs may be directly linked to the mag- 
netic flux of the progenitor star (specifically with reference 
to Ap and Bp stars), in which case the fields would have had 
to survive evolutionary mass loss, if not necessarily Roche- 
lobe overflow. This suggestion is tenuous, and the fate of 
the stellar magnetic field during RLOF is a significant un- 
known in this picture - not just with respect to the initial 
flux, but also the possible operation of a dynamo. We do 



not include any stellar dynamo in the current calculation; 
rather we simply assume that Bg is a constant over time, 
even during mass loss (in §3.2 we discuss the consequences 
of non-constant surface field strengths). 

We use an updated version of Eggleton's s tellar evo- 
lution code (e.g. lEggletonllT97ll: iPols et al.lll995^ to evolve 
systems with initial donor stars of 3, 4, and 5 M© orbiting 
a black hole of mass 7 Mq with an initial orbital period of 
one day. We use the RLOF rate from the previous time-step 
to calculate the magnitude of the irradiation-driven stellar 
wind (eq. 1131 , and from this the magnetospheric radius (eq. 
3J, and finally the angular-momentum loss (eq.jSJ in the cur- 
rent timestep. Overall, this is equivalent to using eq. 1151 . 

Figure |3| shows the outcomes of these evolutionary se- 
quences for a range of stellar magnetic field strengths, and 
with a wind-driving parameter V'/c^ ~ 10~®. With these pa- 
rameters, only systems with donor stars having Bs < 300 G 
will evolve to wide orbits; for larger dipole B fields the sys- 
tems will become compact during the mass-transfer phase. 

In Fig. 13 all the systems begin RLOF with a high mass- 
transfer-rate plateau lasting ~ 1 - 100 Myr (only as long as 
100 Myr for the 300 G stellar fields), at the end of which a 
maximum in the orbital period evolution is reached. Donor 
stars with only 300 G fields reach core hydrogen exhaustion 
after ~ 1 Gyr, and their orbital periods lengthen as they ex- 
pand. The donor stars with stronger fields - having lost mass 
more rapidly - evolve more slowly, continue burning hydro- 
gen in their cores out to > 10 Gyr, and are good candidates 
for short-period black-hole binaries. They have all dropped 
below 1 M0 by ~ 100 Myr after the onset of RLOF; likewise, 
they spend relatively little time above an orbital period of 
~0.5 day s, which agrees well wi t h the orbital parameters in 
table 1 of lLee. Brown fc Wiier3 ll2002l) . At 1 Gyr, with peri- 
ods of ~ 0.25 d and donor masses around 0.5 M©, the IkG 
models produce mass transfer rates of ~ lO"'^" M0yr~^. 

We note explicitly that these systems spend most of 
their evolution with low-mass donors and are most likely to 
be observed as low-mass systems. 

We have also calculated a range of binary sequences 
for less efficient wind-driving energetics, specifically ijj/c? — 
10~^ and 10~*. For these lower wind-driving efficiency fac- 
tors, higher magnetic fields for the donor star are required 
to obtain results similar to those shown in Fig. |21 (see eq. 
18). In particular, Bs must scale as tp'^^"^ . 

These sequences demonstrate that the anomalous mag- 
netic braking scenario is able to produce short-period sys- 
tems for reasonable magnetic field strengths and conserva- 
tive assumptions about the irradiation driven stellar winds. 

3.1 The Long and the Short: Population Statistics 

In order for the proposed anomalous magnetic braking sce- 
nario to be viable for producing short period black-hole 
binaries, it should yield the correct ratio of long to short 
period systems. From the small proportion of A-stars that 
have anomalously high magnetic field s (around 5% - see, e.g. 
lLandstreetlll982i: IShorlin et alj|2002l) . we might expect the 
population of short-period black-hole binaries to be much 
smaller than that of long-period systems, whilst the observed 
populations are roughly equal in size. However, the systems 
driven to lower masses more quickly by magnetic braking 
will live longer than those which are not, and so an imbal- 



6 Justham, Rappaport and Podsiadlowski 



3M0 Donor 4Mg Donor 5Ms Donor 




105 -]q6 -iqS IOIOIOE IQ? -iqS 101[>105 ■] qS -^qJ ■] qS ■] qS ■] qIO 

Years from start of RLOF Years from start of RLOF Years from start of RLOF 

Figure 3. For a nominal wind driving energy conversion efficiency = 10~^ (see text), we sliow the evolution of three different 
companion masses (columns: left to right 3, 4 and 5 Mq) for three different surface magnetic field strengths: 300 G (dashed), 1000 G 
(unbroken line), and 3000 G (dotted). For the lowest field strength, the orbital period finally increases; also higher magnetic fields result 
in higher initial mass-transfer rates, as expected (see eq. 18). The mass-transfer rates have been slightly smoothed for clarity. 



ance in the population birthrates (of around a factor of 20) 
is counteracted by the longer lifetimes of the short-period 
systems. 

A representative lifetime estimate for the short-period 
systems is several Gyr (the 1000 G curves in Fig. I^lall con- 
tinue indicating mass transfer until 10 Gyr, so the lifetime 
depends on choosing a minimum mass-transfer rate cutoff). 
Systems with the same initial orbital period, but with zero 
magnetic field, live only 200-400 Myr (see Fig.QJ. This ratio 
of lifetimes - demonstrated explicitly in Fig.|l]- is approxi- 
mately the required factor of 20 (to compensate for the small 
fraction of Ap stars), and shows that our model is naturally 
consistent with the observed numbers. 

We note that relaxing our simplifying assumption - that 
the magnetic field strength does not decrease as the donor 
is stripped of mass - should make the lifetimes of the short- 
period systems even longer, as in the present calculations the 
low-mass donor stars are still subject to significant magnetic 
braking. More detailed lifetime arguments are also likely to 
depend on how evolved the donors are at contact and hence 
on the post-common-envelope orbital periods. 

Wary of these unknown details, we calculated a range of 
simple population distributions to examine the feasibility of 
matching the observations of these low-mass black-hole bina- 



ries. A detailed presentation of these results is not included 
but we can draw an important general conclusion. The pe- 
riod distribution does seem to be best matched by decreasing 
the surface magnetic field strength when the donor mass falls 
below ~ 1 Mq ; such evolution of the initially strong Ap field 
is eminently reasonable. (We repeat, however, that for all of 
the results shown in Figs. 2-5, the magnetic field strength 
was held constant.) 

We note that - though we expect the stellar surface 
field strength to decay - the onset of transient behaviour 
would naturally reduce the efficiency of this irradiation- 
driven braking mechanism. Equation ^ reveals that the 
torque is proportional to the square root of the wind mass 
loss, and hence (via eq. 113^ to the square root of the mass 
accretion rate. Persistent accretion is thus more efficient 
at driving braking than phases of outburst and quiescence, 
which is an effect that should be considered in more detailed 
future investigations of this mechanism. We also acknowl- 
edge that, in systems where the accretion rate begins to 
exceed the Eddington limit, the irradiation efficiency (and 
hence -i/;) will drop; however, since the mass transfer is it- 
self driven by our braking mechanism, this effect is to some 
extent self-regulating and would not necessarily result in 
widening of the binary. 
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Figure 4. In order to illustrate the contrast in lifetimes for 
magnetic and non-magnetic donor stars, we compare the mass 
evolution of systems with magnetic fields of 1000 G (solid lines, 
as in Fig. 13 and Bs = (broken lines) for initial donor masses of 
3, 4, and 5Mq and an orbital period at contact of 1 d. Nuclear 
evolution and the complete transfer of the stellar envelope termi- 
nate the systems without magnetic braking (which are evolving to 
longer periods) after only a few hundred Myr. The short-period 
systems are still transferring mass after 10 Gyr (for the corre- 
sponding mass transfer rates, see Fig. ISl. 



3.2 Transient Behaviour 

We checked the consistency of our models for compact black- 
hole binaries for their potential behaviour as soft-X-ray tran- 
sients via the familiar thermal-ionizat i on disk instability 
JCannizzo et al.lll98l iKing et alJll996l: Ivan ParadiislllQQd: 
iDubus et al.ll200lt iLasotal 1200 J) At each step in the bi- 
nary evolution we calculated the susceptibihty of the ac- 
cretion disc to instability, either by computing the tem- 
perature that would be attained in the outer portions of 
th e accretion disk due to X-ray irrad i ation (using eq. (8) 
in iRappaport. Podsiadlowski. fc PfahU l2005l) , or by deter- 
mining a c r itical accre tion rate usi n g exp ressions due to 
iKing et alj lll997l) and iDubus et al.1 il999ll . Though these 
different criteria yield considerable freedom of outcome, we 
find that all the methods for determining transient be- 
haviour can match the observations: they agree that systems 
with Porb ^ 0.5 days and Md < IMq are transient (pro- 
vided the magnetic braking at that epoch is not implausibly 
strong) . 

Common to all the criteria for transient behaviour is the 
prediction that during the initial phase of mass transfer - as 
the donor approaches ~ IMq - the systems will be persis- 
tent sources. These systems will be evolving rapidly and thus 
there may be too few of them to have been discovered. How- 
ever, the model predicts that there is a class of black-hole 
binaries with intermediate-mass donors and orbital periods 
< 1 day that are steady X-ray emitters. 

In light of the fact that the wind-driving mechanism 
becomes significantly less efficient in transient systems, we 
point out that at the periods where the systems are pre- 
dicted to become transient the effective temperatures of the 
donors have dropped within the range where 'normal' mag- 
netic braking operates, i.e. the stars are capable of driving 
their own winds. 
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Figure 5. Evolutionary tracks in the donor-mass— period plane 
for donors initially of 4 Mq with initial periods of 18, 24 and 
36 h (dotted lines), as well as for donors initially of 2 Mq with 
initial periods of 24, 30 and 36 h (dash-dotted lines). We also 
indicate the periods of known compact black-hole X-ray binaries 
(arrows on period axis), with crosses marking the earliest spectral 
types given in Casares |200!|). The solid line marks zero-age main 
sequence (ZAMS) donors which would fill their Roche lobes with 
a 7 Mq companion. Note that less evolved donors are closest to 
the ZAMS line, and for the 4 Mq donor more evolved donors 
lead to earlier spectral types, whilst more evolved 2 Mq stars are 
cooler at longer orbital periods but become hotter as the interior 
is further exposed. 



3.3 Observational Test: CNO Abundances 

As well as population arguments, we can use individual 
objects to test the proposed anomalous magnetic brak- 
ing mechanism. Observations that show CNO-processed el- 
ements on the surface of XTE J1118-I-480 (Haswell et al. 
2002; the system has a 4.1 hour orbital period) are strong 
evidence that the donor is the descendent of an intermediate- 
mass star. Under the assumption that only dynamo-driven 
magnetic braking occurs, Haswell et al. (2002) found that 
only a very small parameter space could have produced XTE 
Jl 118-1-480 with such a CNO-processing signature, with a 
donor star of initial mass 1.5 Mq. Our irradiation wind- 
driven magnetic braking model relaxes these constraints 
considerably. 

Moreover, the far UV spectrum taken by Haswell et al. 
appears to lack emission in oxygen as well as carbon. This 
fact is explained most naturally if the progenitor was more 
massive than 1.5 Mq, such that the full CNO tricycle could 
have operated. As PRH suggested, progenitors with such 
large masses would require a new angular-momentum loss 
mechanism - such as the one proposed in this work. 



3.4 Observational Test: Spectral Types 

Figure |S] compares the results of our binary evolution cal- 
culations with unevolved, primordially low-mass donors in 
the Tcff — Porb plane, and marks the periods of the known 
compact black-hole X-ray binary systems. The companion 
stars in these latter bi naries are obse rved to have cool spec- 
tral types (see, e.g. ICasaresj l2005l . where the donors are 
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listed as broadly mid-K type stars; iTorres at al.l ll2004l) dis- 
cuss their spectral type determination for XTE J1118-I-480). 
This forms a strong constraint on intermediate-mass progen- 
itor stars: donors which have been allowed to evolve by the 
onset of mass transfer tend to have effective temperatures 
higher than required at the known orbital periods (see Fig. 

Note that the behaviour of the 4 Mq donor in Fig.|21con- 
trasts strongly with that predicted for the donors in CVs, 
wh ere more evolved donors produ ce later spectral types (see, 
e.g. iBeuermann et alJlTooS iBara ffc & Kolhl l2000l) . 

We stress that the crosses in Fig. |S] mark the hottest 
temperatures consistent with the observed spectral types. 
This difference between theory and observation is hard to 
reconcile even when approximate non-g rey atmospheric cor- 
rections are applied to our models fsee [ Chabrier fc Baraffe 
1997*; 'B araffe et ai]ll998l: iPodsiadlowski. Ifan fc Rappapor' 
2003). At best, we can force the models to approach the tem- 
peratures of primordially low-mass ZAMS stars by having 
the donors begin transferring mass to the black hole early 
on the main sequence. However, even these values of T^s are 
still too high. 

The conversion between spectral type and effective tem- 
perature is non-trivial and it may be possible for our coolest 
models - late-G type donors - to be mistaken for mid-K 
stars, but we consider it unlikely that such a conversion er- 
ror alone could account for the wide discrepancy between 
the bulk of our model s and the observations. However we 
note that iTorres et al.l ll2004l) find that the donor star con- 
tributed only ~ 55% of the light in their observations of 
XTE J1118-I-480 during quiescence and suggest that unam- 
biguous determination of donor spectral types in such sys- 
tems is challenging. 

We did consider one potential mechanism by which the 
secondaries in these systems might appear cooler than our 
models predict. The companion in the black-hole binary 
GRO J 1655-40 (Nova Scorpii 1994, V1033 Sco) is known to 
be polluted with several a-process element s (for a detailed 
investigation, see IPodsiadlowski et al.ll2002fl : such pollution 
would be expected to increase the opacity of the affected 
layers. We investigated both ad-hoc opacity modifications 
for the surface layers and uniform composition changes. Al- 
though we managed to produce closer correspondence with 
the observations, we cannot claim that the improvements 
were very dramatic. 



4 SUMMARY AND CONCLUSIONS 

In this paper we have proposed that the subset of neutron- 
star and black-hole binaries that form with Ap or Bp star 
companions will experience systemic angular-momentum 
losses due to magnetic braking, not otherwise operative with 
intermediate-mass companion stars. We have quantified how 
a magnetically coupled, irradiation-driven stellar wind can 
lead to substantial loss of systemic angular-momentum. We 
have demonstrated with detailed binary stellar evolution 
calculations that the proposed magnetic braking scenario 
involving Ap/Bp donor stars is effective in allowing such 
systems to evolve to short orbital periods (Porb ^ 10 hr). 
In the absence of such magnetic braking, binaries where 
intermediate-mass donor stars transfer mass onto more mas- 



sive black holes inevitably evolve to long orbital periods, i.e., 
Porb weeks. 

Our Ap/Bp magnetic braking scenario has been applied 
to a specific astrophysical problem, namely the formation 
and evolution of compact black-hole binaries with low-mass 
companion stars. A number of previous studies have encoun- 
tered difficulties with explaining the formation of such sys- 
tems (PRH and references within) . This led PRH to consider 
some form of anomalous magnetic braking to prevent sys- 
tems from evolving to long orbital periods. In this work, 
we have presented such a mechanism that requires no new 
physics - and, indeed, is a consequence of previously pub- 
lished mechanisms acting in close binaries. 

One of the conceptual difficulties with the formation 
of short-period black-hole binaries is that a low-mass donor 
star seems incapable of ejecting the common envelope dur- 
ing the prior formation of the black hole (PRH and refer- 
ences within; see also the Appendix). More massive stars, 
including stars of intermediate mass (i.e., 3 — 5 M0), would 
be more likely to be able to eject the common envelope. 
Even under favourable (or optimistic) conditions in which 
a low-mass donor might conceivably be able to eject the 
common-envelope of the black-hole progenitor, stellar de- 
mographics suggest that the formation of binary systems 
with mass ratios of around 20 is very uncommon (see, e.g. 
iGarmanv. Conti fc Massevlll980l) . 

On the other hand, intermediate-mass stars are not 
thought to be subject to magnetic braking. As mentioned 
above, such angular-momentum losses are required to pre- 
vent the system from evolving to longer, rather than shorter, 
orbital periods. We have applied our Ap/Bp magnetic brak- 
ing scenario to this problem in order to solve the difficulties 
with both the common-envelope ejection and the angular- 
momentum loss mechanism. 

We have calculated a sequence of binary evolution mod- 
els that demonstrate the Ap/Bp magnetic braking model is 
successful at reproducing the short orbital periods and low 
donor masses observed for the compact black-hole binaries. 
Moreover, we have shown that, since these systems have 
a longer lifetime than the corresponding systems without 
magnetic braking, the relative populations of compact and 
wide black-hole binaries could be explained in spite of the 
small fraction of A/B stars that are of the Ap/Bp class. Our 
model also helps to explain the evidence for CNO-processed 
material seen at the surface of XTE J1118-I-480. 

As demonstrated in Fig. 5, however, the effective tem- 
peratures of our model donor stars are significantly higher 
than for those of the observed donor stars. This seems to 
be a generic difficulty with any formation scenario that in- 
vokes primordially intermediate- mass donor stars, and is not 
specifically related to our suggested angular-momentum loss 
mechanism. Any new angular-momentum loss mechanism 
proposed to explain the formation of short-period black-hole 
binaries would have to account for the same mismatch in 
spectral types. Hence in developing our model, we may thus 
have established a substantial constraint on the broad set 
of formation scenarios which begin with intermediate-mass 
secondaries. 

The problem results from the fact that even though the 
initially intermediate-mass star loses much of its mass fairly 
rapidly (driven by magnetic braking), the star is still some- 
what evolved chemically (both in He and CNO abundance) 
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by the time its mass has been reduced to ~ 1 Mq. This 
prevents the donor stars from achieving the cooler observed 
values of Tcff • Even in the limit that the intermediate-mass 
star could lose much of its mass within a few of its thermal 
timescales, the subsequent evolution would not automati- 
cally lead to the evolutionary state of the observed donor 
stars, i.e., undermassive and underluminous for a given or- 
bital period. 

Even though the systems that we have generated have 
donor stars with hotter effective temperatures than are ob- 
served, our analysis predicts that such systems should exist 
and could be discovered in the future. The signature would 
be orbital periods between ~ 5 — 15 hours, donor-star masses 
of ~ 0.5-1 Mq, and effective temperatures 5000-7000 K. If 
black-hole binaries with these higher effective temperatures 
are not found, then this would imply that either the donor 
stars are rarely, if ever, of the Ap/Bp variety, that the pro- 
posed magnetic braking mechanism is not as efficient as we 
have calculated, or that the Ap/Bp magnetic fields do not 
persist as mass is lost from the star. 

In a subsequent paper we plan to explore further sce- 
narios for the formation and evolution of compact black- 
hole binaries. These include (i) rapid mass loss of the 
intermediate-mass secondary at the end of the common en- 
velo pe phase; (ii) explosive common -envelope ejection (see, 
e.g. Ilvanova fc Podsiadlowskil l20o3) : (iii) formation of the 
low-mass donors in situ from the remnants of a failed 
common-envelope ejection; and (iv) evolutionary paths of 
low-mass donors that lead to the system properties that are 
observed. 
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APPENDIX: COMMON-ENVELOPE 
ENERGETICS 

Herein we briefly summarise the arguments that imply that 
a low-mass donor star cannot eject the envelope ol a BH- 
producing star. In the standard formulation using simple 
energetic arguments, the expression for the final orbital sep- 
aration at the end of the common-envelope phase, a/, to 
that just at the onset of mass transfer from the primary to 
the secondary, Oi, is, to a good approximation: 

(e.g., Webbink 1985; Dewi & Tauris 2000; Pfahl, Rappaport, 
& Podsiadlowski 2003; Rappaport, Podsiadlowski, & Pfahl 
2005) where Mi, Mc, and Me are the masses of the black- 
hole progenitor, its core, and its envelope, respectively, M2 
is the secondary mass, and tl is the size of the progeni- 
tor's Roche lobe radius in units of the orbital separation. In 
this expression, A is a parameter that gives the binding en- 
ergy, -Bbind, of the envelope of the primary in the expression 
Sbind = GMiA/e/(i?2A), where R2 is the radius of the pri- 
mary. We have assumed here that the efficiency with which 
the gravitational potential energy between the core of the 
primary and the secondary star can be used to eject the 
primary's envelope is unity. For values less than unity, the 
following argument is only strengthened. 

The quantity in parentheses on the right hand side of 
eq. lETJ has a numerical value of 0.022 ±0.002 M~'^ for aU 
primary masses between 25 and 45 Mq which we consider to 
be the main source of the black holes in question. For larger 
primary masses, it only becomes more difficult to remove 
the envelope. We can then write eq. (lAH as: 

a/ ~ 0.011(a,ri)M2A , (A2) 

where the quantity a^rL is just the radius of the donor star at 
the onset of the common envelope phase. For typical progen- 
itor masses of the black holes, we take the maximum stellar 
radius to be ~2300 Rq. Thus, at the end of the common- 
envelope phase the maximum orbital separation is given by: 



The numerical coefficient in eq. (IA4II is good to ±15% over 
a wide range of primary and secondary masses. 

Thus, for initially low- mass secondaries, i.e., M2 < 
1 Mq the A parameter must be > 0.15 even when the com- 
mon envelope occurs at the widest plausible binary sepa- 
ration. For smaller orbital separations at the onset of the 
common-envelope, the values of A required to avoid a merger 
approach unity. Figure 1 of PRH illustrates this binding 
energy parameter for stars capable of forming a BH (see 
also Dewi & Tauris 2000). For the calculations and assump- 
ti ons made by PRH - including win d mass loss according 
to lNieuwenhuiizen fc de JageJ lll990l) - stars of 25 Mq and 
greater have total binding energies such that A < 0.1 for all 
R > 3OO-R0. Note that those values for A include thermal 
terms, and hence assume that all the energy of recombina- 
tion for the ionised species in the envelope can be liberated 
to aid the common-envelope ejection process. 



afc^25(^]( -^1^] A Rq . (A3) 

\MqJ V23ooj?0y " ^ ' 

In order to avoid a merger of the black-hole progen- 
itor core and the secondary at the end of the common- 
envelope phase, the Roche lobe of the secondary must be 
larger than a main-sequence star of the secondary's mass. 
If we use a simple mass-radius relation for low-mass main 
sequence stars, i.e., R2 = 0.85{M2 /Mq)°'^^ Rq, and employ 
the Roche-lobe radius expression for low-mass stars due to 
Paczyriski (1971), we find the following lower limit for the 
parameter A in order to avoid a merger of the core and sec- 
ondary: 

A>0.15r^^)^ . (A4) 

\ ^Imax / \/M2 



